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Abstract: We consider a ‘color density matrix’ in gauge theory. We argue that it systemat¬ 
ically resums large logarithms originating from wide-angle soft radiation, sometimes referred 
to as non-global logarithms, to all logarithmic orders. We calculate its anomalous dimension 
at leading- and next-to-leading order. Combined with a conformal transformation known to 
relate this problem to shockwave scattering in the Regge limit, this is used to rederive the 
next-to-leading order Balitsky-Fadin-Kuraev-Lipatov equation (including its nonlinear gen¬ 
eralization, the so-called Balitsky-JIMWLK equation), finding perfect agreement with the 
literature. Exponentiation of divergences to all logarithmic orders is demonstrated. The 
possibility of obtaining the evolution equation (and BFKL) to three-loop is discussed. 
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1 Introduction 


Collimated sprays of particles, or jets, figure prominently in high-energy collider physics. 
This has led to a growing interest in the characterization of jet shapes and event shapes, 
with the goal to extract as much information as possible about underlying hard scattering 
events. The pencil-like nature of jets implies that one often encounters disparate angular and 
energy scales. These lead to large logarithms in theoretical calculations, whose resummation 
is necessary to obtain controlled, precise predictions. Theoretically, in analytic studies these 
large logarithms are often the only terms which one may hope to predict in an amplitude or 
cross section at higher orders in perturbation theory, and thus could potentially help reveal 
new structures. Both of these reasons make them especially important. 

Thanks to developments spanning many years, resummation for most observables of 
interest is now possible. In the case of so-called global observables, which involve complete 
(‘global’) integrals over final state phase spaces, a critical ingredient the exponentiation of 
infrared and collinear divergences [1-6]. This predicts in a quantitative way the logarithms 
left after the cancelation of infrared and collinear divergences, guaranteed on general grounds 
by the Kinoshita-Lee-Neuenberg (KLN) theorem [7, 8]. There exists however non-global 
observables, for which phase space cuts lead to soft radiation not being integrated over all 
angles (‘not globally’), which are considerably more difficult to resum [9, 10]. 

The aim of this paper is to propose a comprehensive theory of non-global logarithms, to 
all logarithmic orders and finite Nc- This theory will be closely related to the BFKL theory, 
which controls another class of difficult logarithms in the Regge limit (high-energy scattering 
at fixed momentum transfer) [11, 12]. 

To set the stage we consider a generic weighted cross-section of the form 

0 - = ^ / dUn\AQ^n{pi,---,Pn)\‘^u{{pi}) (1.1) 

n 

where dUn is the phases space measure for n partons. The measurement function u{{pi}) 
specifies the details of the measurement, including various vetoes etc. For suitable infrared- 
and collinear- safe choices the cross-section will be finite order by order in perturbation theory. 
As a preliminary simplification (to avoid initial state radiation), in this paper we will assume 
that the initial state is a color-singlet state of mass Q, and assume massless final states. 

A time-tested strategy to resum large logarithms is to introduce intermediate, divergent, 
matrix elements, which are to be renormalized at a suitable factorization scale. The template 
is Wilson’s operator product expansion, which expresses the short-distance limit of a corre¬ 
lator in terms of short-distance OPE coefficients, anomalous dimensions, and long-distance 
matrix elements. The factorization scale p, whose dependence is controlled by anomalous 
dimensions, cancels between the OPE coefficients and matrix elements and provides a han¬ 
dle on large logarithms. Our main proposal is that the pertinent operator for resumming 
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non-global logarithms is the color density matrix: 


<rlU] ^ Y. 


dUr 




We call it a density matrix because it is linear in both the amplitude and its complex conjugate. 
It is a functional of a continuous field of unitary matrices which depend on a two- 

dimensional angle and live in the adjoint representation of the gauge group. Pictorially, U, 
shown in fig. 1, is a (local) color rotation between the matrix element and its conjugate. A 
closely related construction has been used to describe parton showers at finite W [13]. 

The physical motivation for eq. (1.2) is that the information carried by cr[U] is clearly 
necessary to fully characterize the distribution soft gluons. Since soft radiation can be trig¬ 
gered by any other colored parton with a higher energy, keeping track of color charges in every 
direction, as cr[U] does, indeed seems unavoidable. The information in cr[U] is also intuitively 
sufficient: due to coherence effects, soft gluons are affected by the color charge carried by 
harder partons but generally not by other details. 

Contrary to the original weighted cross-section, the density matrix cr[U] is infrared di¬ 
vergent. We propose, and will demonstrate, that these infrared divergences exponentiate in 
terms of a well-defined anomalous dimension. This supports our claim that the information in 
is sufficient. After cancelling these divergences (see eq. (2.7)), the ‘renormalized’ density 
matrix then depends on a factorization scale as 

p] = K{U, 6/6U, e) p] . (1.3) 

The anomalous dimension or Hamiltonian K assumes the form of a functional differential 
operator. Its one-loop expression, computed in eq. (2.14), repdoduces earlier formulas derived 
in the literature to deal with non-global logarithms [10, 14]. 


A rA. 

^ dp das 


1.1 Structure of the resummation 

For concreteness, we briefly sketch how the formalism would be applied in a complete calcula¬ 
tion of a specific non-global event shape, the (cumulative) hemisphere mass function, at finite 
Nc- The hemisphere mass function is the distribution of the invariant masses in the left- and 
right- hemispheres of a decaying color-singlet state of mass Q (for example, a Z-boson). In 
the di-jet limit mL,mR <C Q, the observable simplifies to 

n{x,Y,Q) = ^ IdUr^lAg^^fe |^x-^0i(f)(pO+pf)^0|^y-^0K(f)(pO_pf) 

(1.4) 

Here OiiJi{i) = 0(=Fpf) project onto respective hemispheres. Logarithms of X/Q and Y/Q are 
global and can be resummed using standard techniques (see ref. [15] and references therein). 
When A <C T <C <5, 77 receives additional logarithms log(T/A) due to the phase space cut 
between the two hemispheres. These are the non-global logarithms. 



- 3 - 










Figure 1. Color density matrix. For each colored final state, an independent color rotation is applied 
between the amplitude and its complex conjugate. 


To resum them we take Y as the hard scale and X as the soft scale. At the hard scale 
we drop the X-dependent step function, and define a density matrix involving only Y: 


a[U;Q,Y] = Y, / dn„ 




(1.5) 

This is of the form of eq. (1.2). Infrared divergences introduced by the U matrices exponen¬ 
tiate and renormalizing at the scale p = Y (see eq. (2.7)) removes large infrared logarithms. 
Concretely, in perturbation theory, is polynomial in f/’s and can be viewed as a book¬ 
keeping device encoding the orientations of radiated particles. It starts with the dijet term 


a^^'^[U;Q,Y,p = Y]=TT[U{n)UHn)] + 0 (^ 2 ), ( 1 . 6 ) 


with n, h null vectors along the positive and negative z directions. To resum logarithms 
log(y/X) we now use the RG equation (1.3) to run down to the scale X, where we deal 
with the infrared part of the measurement. In the leading-log approximation, the 6{X — ...) 
factor in the observable effectively removes left-hemisphere radiation between the scales X 
and Y (except within a cone of small radius R ~ -sJXjQ centered around the left jet). The 
IR measurement can thus be phrased in terms of an averaging 


7^(X, y, Q) = (^’-^“[17; Q, y, M = ^])iR (1-7) 

where ^ 

+ 0(5^) with (t/“''(0i))o = <5“%(i). (1.8) 

The step function Q'r{i) = 0i?(i) + 6{—Pi/p^ + cosR) allows real radiation inside the right 
hemisphere and left cone. It is important that the averaging depends only on angles, since at 
this stage information about energies is lost. 

^ Even though the right-hand-side is not a unitary matrix, the equation makes sense since the average of 
unitary matrices need not be unitary. 
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In principle both the UV and IR measurements (1.5), (1.8) are to be computed order by 
order in perturbation theory. Resummation is considered successful if theses series contain no 
large logarithms. Since the Hamiltonian K is only meant to resum soft logarithms between 
X and Y, the UV and IR measurement functions actually do contain further logarithms of 
Q and R, but these should be resummed using other, conventional, means. Comparison of 
this procedure with the leading-log prescriptions of refs. [9, 10, 16] is discussed in section 2. 
In keeping with the logic of factorization, in this paper we will concentrate on the (universal) 
evolution K and leave calculations of (process-dependent) measurement functions to future 
work. 

A remarkable fact about K is that it is also (essentially) the Balitsky-Fadin-Kuraev- 
Lipatov (BFKL) Hamiltonian, that is, the boost operator of the theory in the high-energy 
limit. The same Hamiltonian K thus simultaneously governs non-global logarithms and the 
Regge limit. This was observed mathematically from the one-loop expressions in refs. [10, 
14]. A general explanation is available, however, using a conformal transformation, which 
extends to higher loop orders [17]. One thus anticipate the difference in QCD to be at most 
proportional to the /3-function. 

Since this correspondence will play an prominent role it is helpful to include a rough 
explanation here. High-energy forward scattering (for example the elastic pp —)• pp amplitude) 
amounts to taking an instantaneous snapshot of a hadron’s wavefunction, so pictorially it 
measures the amplitude for a virtual shower to form inside the hadron and then recombine. 
This is illustrated in fig. 2(a). This is also roughly what the density matrix ) [U] of a 

decaying virtual hadron measures. Importantly, however, one measurement is instantaneous 
while the other takes place at infinity. To relate them requires the nontrivial conformal 
transformation of ref. [17]. 

In this correspondence with the Regge limit, the color rotations in (j\U] implement the 
shockwave of the Balitsky-JIMWLK framework [18-21]. Here the ‘shockwave’ is inserted 
at infinity between the matrix element and its complex conjugate. This was our original 
motivation for defining eq. (1.2). (Mathematically similar considerations were used in refs. [22, 
23] to exploit the conformal symmetry of the BFKL equation.) The UV and IR measurements 
represent, respectively, the projectile and target impact factors. 

The aim of this paper is to analyze the properties of the Hamiltonian K and to calculate 
it explicitly to the next-to-leading order. The lessons learned from this calculation will then 
lead to an immediate proof of exponentiation. As a cross-check of the calculation we will 
compare against results obtained in the context of Regge limit scattering. 

This paper is organized as follows. In section 2 we review known facts regarding the 
exponentiation of infrared divergences and factorization of soft emissions. We illustrate the 
formulas by giving the leading terms in perturbation theory for the various ingredients. We 
also verify that the procedure around eq. (1.8) reproduces the established resummation of 
non-global logarithms at leading-logarithm order. In section 3 we perform the two-loop 
calculation. A key finding will be the possibility to express all terms in K as finite integrals 
over well-defined, finite and gauge-invariant building blocks. The final result is recorded 
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Figure 2. (a) Scattering in the Regge limit. The thin shock is the Lorentz-contracted target, (b) 
Branching of soft gluons. To connect the pictures one ‘folds’ along the target and sends it to infinity. 


in subsection 3.6. In section 4 we compare our result for K against the two-loop BFKL 
equation. We will find perfect agreement in conformal theories, with, as expected, a relatively 
compact correction term proportional to the /3-function in QCD. In section 5, using the lessons 
learned from the two-loop calculation we derive formal expressions for K at three-loop and 
beyond, and demonstrate exponentiation in general. Conclusions are in section 6. A technical 
appendix reports complete details of our evaluation of the real-virtual contributions at two- 
loop. 


2 Conventions and review 


To set our conventions we now review the exponentiation of infrared divergences and give 
explicit formulas for the relevant objects at one loop. We also discuss the resummation of 
non-global logarithms at leading-log. 


2.1 Soft factorization 


As described in refs. [4-6], the exponentiation of infrared and collinear divergences is con¬ 
trolled by a soft anomalous dimension: 


An = V exp 



(A, q:s(A), e) 


X as(Ai),e) • 


( 2 . 1 ) 


For a gentle introduction we refer to ref. [24]. The infrared-renormalized amplitude Hn is 
also called the hard function and is finite (in this paper we use only ultraviolet-renormalized 
amplitudes). It depends on the choice of a factorization scale //, 


3^ 


dfj, 


+ /5(q^s) 


d 

das 


Hn{n,as{n),e) 


Tnin, as{^^), as(/^), e) • 


( 2 . 2 ) 


It is important to note that since F„ acts as a matrix in the space of color structures, the 
path-ordering symbol cannot be omitted. 
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We work in D = 4—2e dimensions and the coupling constant depends on scale through 


d 

= I3{as,e), 


/3(as,e) = -2eas + /3(as) 


(2.3) 


At one-loop /3{as) = —2^bo with bo = AtifTf nsTs ^ theory with np flavors of Dirac 
fermions and ns complex scalars (in QCD Ca = 3 and Tp = j)- The solution is 


-2e r 


asin) = as(^o)^^ 


1 + 


boasino) I - fj. ^ 


27r 


2e 


(2.4) 


The integral in (2.1) converges and produces the desired 1/e poles provided that e is negative 
enough that the coupling vanishes in the infrared. 

In the literature T^ is often written as being e-independent, which defines minimal sub¬ 
traction schemes. We keep the more general notation since it will be useful to consider 
non-minimal schemes. As long as T^ remain finite as e —)• 0, different schemes are related by 
finite renormalizations of Hn. 

Amplitudes with m soft gluons are also known to factorize in a simple way [6, 25], 

(Ai, as{^l),e) as(M), e) x a^iiJ.), e), (2.5) 

up to powers of k provided that all of {ki,..km} are softer than other scales in M. We 
have only indicated the color and Lorentz indices of the soft gluons (to be contracted with 
polarization vectors e/'). Since is the same as in eq. (2.1), this formula states that soft 
gluons can be ‘tacked onto’ an amplitude without recomputing it. 

Similarly to T^, the soft currents Sm are matrices in the space of color structures. Ac¬ 
cording to eq. (2.2) they are finite and have factorization scale dependence 


d d 


-SmifJ-, Osifi), e)r„(Ai, as{n),e). ( 2 . 6 ) 

Our main proposal is that the color density matrix admits a similar factorization. 


a[U] = V exp 


d\ 

-A(A,as(A)) 


fT’'-[[/;/r,«s(/r),e] 


(2.7) 


where is finite and obeys the RG equation (1.3) quoted in the introduction. Fur¬ 

thermore K is independent of the measurement function u{{pi}) appearing in the definition 
of (t[U]. These will be shown in section 5 to be essentially combinatorial consequences of 
eqs. (2.1), (2.5). 

2.2 Leading-order expressions 

The tree-level emission of one soft gluon is controlled by Weinberg’s well-known soft current: 


= 9'^ RiSi (fci) + 0{g^) where S'/ {k 


I = 


A 


(di-ki ’ 


( 2 . 8 ) 
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Figure 3. One-loop evolution of the density matrix, (a) Real emission of one soft gluon. This adds 
one [/-matrix (cf. eq. (2.11)). (b) Virtual correction. 


where is the operator which inserts a color generator on leg i and = (l,d))^ is a null 
vector proportional to p*. These obey i?^] = and our normalizations are such 

that Tr[r“T^] = and Tr[l] = Nc- The soft anomalous dimension at one-loop is [4] 

= -2 E + E (2.9) 


where the collinear anomalous dimensions are = —6o for gluons and = —SCp for 
quarks. We loop-expand using the uniform notation: r„ = with 


r(l + e)r(l-e) 
r(l - 2 e)( 47 r)-^ ' 


To find K at one loop it suffices to compute (j\U] to that accuracy and compare the 
divergence with eq. (2.7). The real emission contribution to <j\U] has an infrared divergence 
when an additional gluon is emitted at a wide angle, as shown in fig. 3(a). It is given as the 
square of the soft current (2.8): 


'‘■“’teal = E 




dvr 




where is the (infrared) pole from the energy integral, 

/■M 




daai-^^5f(a/3o)5f(a/3o) 


a 




i/[2e] aoiaoj 


( 2 . 11 ) 


( 2 . 12 ) 


with aij = ^2 Here the sums run over the [/-matrices present in (j\U] (which 

at finite order in perturbation theory is a polynomial) and we use the abbreviation = 
The notation here is slightly different, but compatible with, that used in eq. (2.8): 
is the operator which replaces Ui with UiT°‘. Similarly L“, representing the color charge 
in the complex conjugate amplitude, replaces Ui with T^-Ui. These obey: 


[R';,Rl] = ir%kR-, 




[R, L]=0. 

















The virtual corrections (fig. 3(b)) generate products of the type LL and RR with no 
extra U. An important constraint is that = 6 °'^] must be evolution-invariant, since 

this correspond to the total cross-section which is finite by the KLN theorem. Equivalently 
K must vanish when U is the identity field. This unambiguously determines the LL and RR 
terms. Using the identities 


jjaa 
^ i 


rya _ j a 




(2.13) 


which in particular yield Lj = i?j when = 5“^, the (unique) solution is easily seen to be^ 




^(1) 

dvr 


C7o““ {L'}R'^ +R^ 


L a\ T a T a -na j^a 

j) 




(2.14) 


This gives the complete scale dependence of the density matrix including non-planar 

effects (and therefore, by the expected factorization, any non-global logarithm at leading-log). 

We review a few known facts about this equation. 

• Taking the’t Hooft planar limit Nc ^ oo with A = g'^Nc fixed, eq. (2.14) becomes for 
the dipole Uij = ■^Tr[UiUj]: 

KUi2 = ^ j (2 CfUi 2 - ^Tv[T-UiT-'ul]ur'\ + 0 (/) • (2.15) 

lOvr^ J dvr aoiao 2 \ Ac J 

Using simple color identities this reduces to a closed nonlinear equation:^ 

KUu = J^[ (Ui 2 - U 10 U 02 ) + 0(A2, 1/W) . (2.16) 

OTT^ J dvr aoi «02 

This is the Banfi-Marchesini-Smye (BMS) equation governing non-global logarithms 
in the planar limit [10].^ Let us be more precise. As stated in the introduction, the 
functional RG equation (1.3) is to be integrated from the UV to the IR, starting from 
e.g. the dijet initial condition cT[f7] = Unn (1-6). In the IR one performs the average 
(1.8). In the planar limit, the averaging reduces to evaluating the functional at one 
point, cr\Uij = d):j(i)0):j(j)], so the procedure is equivalent to evolving the argument of 
the functional from the IR to the UV, e.g. the function Uij. This is how the relevance 
of eq. (2.16) arises. 

^ A term ■ fij (!/“!/“ — Rt Rj) would also satisfy the KLN theorem and preserve the reality of a provided 
that its coefficient is imaginary. The imaginary part of the explicit expression (2.9) however shows that fij oc in 
is constant and thus cancels out using color conservation. 

^We have used: = U^T‘^Uo and Tt[T‘^XT‘^Y] = |Tr[A]Tr[y] - ^Tr[Xy]. 

^ In addition, compared with ref. [26] which deals with the hemisphere function, one needs to set Ui 2 ^^ = 
6)j(l)@jj(2)[/i2®'^®; the step-functions factors are stable under evolution. At leading-log collinear divergences 
exponentiate independently so the R term in O'n in eq. (1.8) does not interfere with the non-global part. 
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• Away from the planar limit, eq. (2.14) coincides with the generalization of the BMS 
formula derived in ref. [14]. Again, as in the BMS case, the two forms differ only by 
multiplication of the U matrices by step functions, which commute with the evolution. 
The averaging procedure (1.8), done in the infrared, is as in refs. [14]. 

• The work of ref. [14] was at least partially motivated by analogy with equations de¬ 

scribing the Regge limit. Using the substitution given below in eq. (4.2) (namely, 
aij —>■ (xi—Xj)‘j_ and f —)> ^-^), eq. (2.14) is indeed recognized as the Balitsky- 

JIMWLK equation [18-21]. 

• The double-sum notation in eq. (2.14) is most natural in a perturbative context where 
cr[[/] is a polynomial in the U’s. Since the evolution increases the number of t/’s, for 
solution it can be better to view AT as a functional differential operator acting on cr[U]. 
This is achieved by the following simple substitutions (done after normal-ordering all 
L, R’s to the right of U’s) [27, 28]: 

( 2 . 17 ) 

These and obey the same commutation relations as those defined previously, and 
in fact after substitution into eq. (2.14) one finds the same action on any polynomial 
cr[[/]. This reveals eq. (2.14) as a functional second-order differential equation of the 
Fokker-Planck type, whose solution can be importance-sampled via lattice Monte-Carlo 
techniques [27, 28]. For studies of 1/A^c effects in the context of non-global logarithms, 
see [16, 29]. 


• Also well-studied is the weak-field regime where all matrices are close to the identity. 
Following ref. [30] and references therein, we write Uj = and expand (2.17) in 

powers of W. This can be streamlined using the Baker-Campbell-Hausdorff formula, 
which gives 


igL^,,igR°: = 


5Wf 


±-/ 

2-’ 


abx 


_ 


_)_ ^ Jeby^yx^ry 


12 


J 5W^ 


720 ^ ’ 


(2.18) 


Plugging this into the one-loop Hamiltonian yields after a small bit of algebra (only the 
first two terms contribute) 


kW 


ebb c 


d?Q.Q d?Q.i(PdQ.2 ai2 
dvr aoi«02 


(VFf'-lTo“')(^l' 



52 

5Wf5W!^ 


(2.19) 


up to nonlinear terms of the form 6 K ~ g^W'^ . This is one form of the one-loop 

BFKL equation and its (‘BJKP’) multi-Reggeon generalization [11, 12], valid for color- 
singlet states (see ref. [30] and references therein). It acts on functionals a[W] where 
is identified as the Reggeized gluon field. This identification will play a useful role 
later in this paper. 
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(«) 



(b) 


Figure 4. Building block for next-to-leading order computation: amplitude for two soft particles. 
Solid lines are eikonal Wilson lines, (a) Two soft gluons. The non-abelian part of the first graph gives 
a connected contribution, (b) Two soft fermions or scalars. 


• Finally, we did not prove in this subsection that divergences do exponentiate according 
to eq. (2.7). We simply read off the exponent from a one-loop fixed-order calculation. 
Proofs to leading-logarithm accuracy are in refs. [10, 14] and an all-order demonstration 
is given in section 5. 


3 Evolution equation to next-to-leading order 

We now present a calculation of K to the next-to-leading order, by matching two-loop infrared 
divergences in cr[U] against eq. (2.7). The computation will be phrased exclusively in terms 
of convergent integrals over building blocks with a clear physical interpretation (renormalized 
soft currents), which will shed light on the exponentiation mechanism. We perform the 
computation in a general gauge theory, although at intermediate steps we only write formulas 
for color-adjoint matter. The reader not interested in the technical details can skip directly 
to the final result in subsection 3.6. 


3.1 Building blocks: soft currents 

A natural building block is the tree-level amplitude for emitting two soft gluons. It can be 
written naturally as a sum of disconnected and connected contributions: 


Sl^{h)S^{k2) +g^Yl ir’'''RlSr{h,k2) + 0(/), (3.1) 




with S^{ki) 


the one-gluon soft current introduced previously. The connected part 


Sriki,k2) 


1 

2(3i-{ki+k2) 


_l3i-ki 


A-^2 krk2 


(3.2) 


follows directly from the Feynman graphs shown in fig. 4(a) [25]. To optimize the notation 
all color generators are implicitly symmetrized: RfBlj —?• which is relevant when 

i = j. This notational convention (borrowed from ref. [31]) is why the connected part is 
proportional to 
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Figure 5. Second building block: one-loop soft current. 


To familiarize ourselves with the notation we review the transverseness check: = 0. 

For the individual Si one finds 

k,,snh) = 1. h,sr(h,k) 4 (^ - |h) . (3.3) 

We need to use color conservation in the form of the identity Yli — 0- Since this holds 
when Ri is inserted to the right of an operator product, the implicit symmetrization in 
eq. (3.1) produces commutators. For example the divergence of the first sum is 

^ R-R) Sf (/ci)5J(fc2) = \{Rl R]} S-{k2) = 

i,j id 


Y. «?isr(fe) = E . 

i i 


This is easily seen to cancel the second term in the parenthesis in (3.3), up to a /Sj-independent 
term which itself cancels due to Yi — Oj proving transverseness. 

Pairs of soft fermions or soft scalars can also be emitted (fig. 4(b)). For notational 
simplicity we carry out all intermediate steps in a theory with color-adjoint Weyl 

fermions and real adjoint scalars (the hnal result will be trivial to generalize). Then: 


Sff{k,, k2) = ig^ A:2), S^^iki, k^) = ig^ Y. ^ 2 ) (3.4) 


with 


s(f{kiM) = 




2Pr{ki+k2)kvk2' 2^iiki+k2)krk2' 

The second natural building block is the next-to-leading order soft gluon amplitude 
We take it from the literature and quote the result in (A.l); representative graphs are shown 
in fig. 5. 

Equation (2.7) gives the next-to-leading order kernel as the divergent part of the following 
combination: 


7^{2)^(0)ren 



kW 

2e 


^(l)ren 


1 

2 



+ 


boK(^^ 

4g2 



1/N 


We will now see that this can be expressed in terms of soft currents. 


(3.6) 
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3.2 Double-real emission 


We begin with the terms in the NLO kernel which involve two wide-angle soft partons, and 
thus generate two additional U factors. The double-real contribution to is by definition 
(suppressing color indices) 


a 


( 2 ) 


double real 


7r(27r)“2ecp 7r(27r)“2ecp 


UoUo> 


|^„+2|^(o/3o, b/3o') • 


0<a<6<oo 


(3.7) 

(The integrals have compact support due to the momentum-conserving J-function in A, and 
we do not show a factor dllnu{{pi}) associated with the underlying hard event.) The trick 
to evaluate (3.6) is to find compatible integral representations for and For 

we already have eq. (2.12) and subtracting it from leaves simply 


(l)ren| 

heal 




7r(27r 


i-2e 


cr 


-Uo 


^''da\An+i\ ia/3o). 


(3.8) 


The essential point here is that the matrix element factorizes in the soft region, |j4„_|_ip(a/3o) —>■ 
\SnAn\'^, SO that subtracting is equivalent to removing the integration region a < p 

(to all orders in e). Invoking factorization similarly, eq. (3.6) can be re-written as: 


0<a<6<oo 0<a<fi 0<a<b<fi 

{i<b<oo 

= I (^F{a,b)-Fia,b)l^^^)+ j F{a,b). (3.9) 

0<a</i /j.<a<6<oo 

a<b<oo 


Here F{a,b) denotes the integrand in (3.7). This formula is also exact in e. One can see 
that the second integral is finite and the first integral has no subdivergences. (Except from 
collinear regions, which are dealt with in the next subsection.) After scaling a —)• aft in the 
first integral to extract the pole, one thus just get: 


i^(2)| 


double real 


- UoUq' I ada 

TT TT Jo 


(s{a/3o, l3o/)S{af3o,(3o^)* 



(3.10) 


This is the desired formula, which expresses the double-real contribution to as a conver¬ 
gent integral over tree-level soft currents. The integrand measures the extent to which two 
soft emissions are not independent from each other. 

Using the explicit expressions (3.1) the formula yields two nontrivial color structures 


X(2)| 


double real 


E 

i,j,k ' 




dvr d-TT 


+E / ^ (i?<+ Rt'^i) (3.11) 
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(shown in fig. 6). These multiply angular functions 


^S;00'=-8 J ada (5f (a/3o)S’)'(/3o0^r(«/5o,/3o') + S^{/3o)S^{aMSi:''i/3o,aM - 

^ij- 00 '=~^ / ('S'y (a/3o, /3o/)5j‘‘'(a/3o, M + a/3o^)S^’'{l3o, a/3o/) + matter - ■ 

We note the absence of a fully disconnected (Abelian) color structure: since its squared 
amplitude is proportional to l/o^ it disappeared before integration. The matter contribution 
in the last parenthesis is, in full: + S(^. The a-integrals are 

elementary (and convergent!) and after some straightforward algebra starting from (3.2) yield 


AT, 


( 2 ) _ 


1 


aoiOo' 


J L 


a 


i] OiikO^jk 


otojO-ik ajkao'i 


OiQO' 


^(2) 


a 




aoiaoo'ao'j 

+ 


1 + 


a0fc«0'fc OlQkOlOO' 
OOO'Wj 


aoo'«o'fc 


log 


a, 


Ok 


a, 


(3.12a) 


ao'jCioj 


O'k 

adj 


iQg (^QjO^O'j ^ (^Weyl 4) “b log 


aoo' aoi«o'i “ «o'Woi 


(nr - 2 n. 


adj I 
Weyl 


2q!qo/ 


+ ao'iOoj aoi«o'i 2 

aoiao'j — ao'j^oj ao'i«Oi 


(3.12b) 


/o') 

For A'b.QQ, we have used symmetry in (i -H- j) to simplify. The expression is especially 
compact in AA = 4 SYM (the first term). The second and third terms represent, respectively, 
an (adjoint) M = 1 chiral multiplet and a scalar. The rational structures are such that 
all potential divergences associated with the /3i,/3j and /3fc regions cancel. There remains a 
divergence as / 3 o —>■ /3q', proportional to the gluon collinear anomalous dimension 7^ oc bo, 
which will be canceled shortly. 

At this point we could stop: using simple and not-so-simple physical considerations one 
could determine the full result using only what we have so far. We find it instructive, however, 
to continue with the explicit computation. 


3.3 Single-real emission 

We now turn to terms with one radiated parton; since these contain only one U field at a 
wide angle, these will combine with and cancel the collinear divergences in eq. (3.12). 

A good starting point is to express the single-real contribution in terms of the hard 
function (2.1): 


•^single real — 


d‘^-^^Qo 

47r(27r)“2'^cr 


daa 


l-2e 


Ve-fof^-+^Hn+i{al3o;fi) 


(3.13) 


In the soft region we can replace the amplitude by a soft current times Hn. It is useful to run 
the soft current to its natural scale, ^ = a, using eq. (2.6). The energy integral then becomes. 
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Figure 6. Sample two-loop diagrams included in eqs. 


(3.11) and (A.7). 


formally to all orders in perturbation theory, 
/■M 


daa 


l-2e 


■pe -^0 ^ ^'‘+^Si(aido; a) 

roo 

+ / daa}~‘^^ Ve~-f^o Hn+iia^o) 

J U 


TM dX T 

Ve-J- 


(3.14) 


Comparing against the factorization formula (2.7), since r„ ~ -^Ivirtuai that the second 

integral represents a finite correction to the finite coefficient The a-integrand on the 

first line, on the other hand, is nicely identified as the following shift to the exponent: 


K 


d^nr 


single real 


TT 


Uo\Si{n/3o; fj.)[ 


(3.15) 


This gives the single-real contribution to K, formally to all loop orders. This generalizes in 
the simplest conceivable way the leading-order result: one simply evaluates the loop-corrected 
soft current with energy set equal to the renormalization scale The ‘bar’ on S accounts 
for a collinear subtraction, to be defined shortly (as well as for an 0(e) mismatch between 
r„ and I virtualI^6cause the soft current is used at its natural scale, the series for (3.15) 
contains no large logarithms and the bo term in eq. (3.6) is automatically accounted for. 

We now deal with collinear divergences. We do so by subtracting suitable splitting 
functions. First, in the preceding subsection, we replace eq. (3.8) by 


/ j 2—2eQ roo ( 

47r(27r)-2°cr J \^^o\An+i{a/3o)\ - ^6»(fe-a)|Spliti(a/3o; fei)| |Au| 

- a) -h ^6l(/i-a)|Splitj((fci - a)bo; ki)\‘^\An\‘^\^^A ■ 

(3.16) 

The next-to-last term removes the a —)• 0 limit of everything to its left. The formula differs 
from eq. (3.8) only by the splitting functions, which are added in such a way as to introduce 
no new soft divergences. The function Splits(pi; ki) (representing the amplitude for particle i 
to split into two with momenta a/3o and ki — o/3o, symmetrized between the two), is required 
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to have the same integrand-level collinear singularity as when /3o || ki, ensuring 

convergence. This is guaranteed to exist by the factorization of amplitudes in collinear limits, 
and an explicit expression is given in (A.3). 

We stress that the subtraction is not written as an integral over the phase space of 
(n + 1) particles with the original total momentum, since a and /3o do not enter An- This 
would probably be very inconvenient for numerical integration, which is a limitation compared 
to more general methods, for example, dipole subtraction [32], However this is simpler and 
will suffice for our purposes. 

To avoid changing one needs to make a compensating subtraction in the virtual 

contribution: ~ 2 'Re[Hn '^where 


fjA) = 4 ( 1 ) _ 

-IJ. tf. 


da a 


-2e 




Iq a h J 47r(27r) ^'^cr aoiaoj 




4 ( 0 ) 

-^-rj 




da a 


l-2e 




-2e 


47r(27r)“^'^cr 


inCi 


|Spliti(a/3o; ki)\ A® + ^. 


(3.17) 

This differs from the naive subtraction of the virtual part of (the R^Rj term), which is 
essentially T^ up to 0(e) terms, by terms involving the splitting amplitude and compensating 
those in eq. (3.16). The zvr term is such that, as verified in appendix, Hn'^ — Hn'^ is finite. 
The barred hard functions can thus viewed as simply the hard functions in a different scheme. 

What is important is that the same subtractions can be applied to all expressions in the 
preceding subsection. For example, in the double-real term in eq. (3.7) one simply subtracts 
(|Splitg(a/3o, 6/3o')|^+Z]i |Spliti(a/3o, fcj)|^) |^n-i-i(&/So')|^ from the integrand. In the T^^) terms 
one subtracts only the soft limit. In this way all two-particle-collinear divergences are removed 
from the preceding subsection, at the only cost of adding a piece to eq. (3.10): 



f dQp dflp' ada 

double real — / ^0 1 ^ , 

linear in 17 J VT TT J q \ a 

Q 

Sphto(o/3o', {l + a)l3p)\^ 

5(o)(/3o)|'- 

l/a2^ 

(3.18) 


This makes it collinear-safe when added to the integrand. The argument of the splitting 
function is such that most energetic particle of the pair has energy /r (which we scaled out). 
Regarding the virtual corrections, the only effect of (3.17) is to put the bar on S in (3.15). 

Both eqs. (3.15), (3.18) are explicitly finite and to be evaluated directly in four dimensions. 
Their evaluation is conceptually straightforward and detailed in appendix A. It turns out that 
the result could have been anticipated using (not so trivial) physical considerations, so here 
we concentrate on explaining these considerations. 

The most critical consideration is gluon Reggeization, or, more broadly, the connection 
with BFKL. As mentioned in introduction, K is the BFKL Hamiltonian in disguise (up to (3- 
function terms). Interactions between Reggeized gluons are constrained by physical principles 
such as Hermiticity of the boost operator and signature conservation U°'^ —>■ which are 

not self-evident from the perspective of non-global logarithms. 
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We consider the weak-field regime where U = with the ‘Reggeized gluon’ field 

W small. Linearizing the Hamiltonian yields Reggeon-number conserving terms given in 
eq. (2.19), as well as 2 —)• 4 transitions (between states with with different powers of W) 
at order g^, and so on. Hermiticity of the boost operator (with respect to the specific inner 
product given by the scattering amplitude of left- and right- Wilson lines) then predicts 4 —?• 2 
transitions at the same order, whose existence is indeed well-known [33]. These are the terms 
which close the so-called Pomeron loop. Now when reverting to the current power-counting 
which treats {U — 1) as 0(1), instead of 0{g), these 4 —)• 2 transition become a three-loop 
effect (see ref. [30] and references therein). Signature forbids 3—7-2 transitions. Hence the 
remarkable statement that K must be triangular at one- and two-loop [30]: 


(fnoW^ 






>0 (L = l,2). 


(3.19) 


Mathematically, this formula (just with the L = 1 case) can be seen as equivalent to gluon 
Reggeization, since it ensures that sectors with different powers of VP’s can be diagonalized 
independently at one loop. One then expects the Reggeized gluon (VP held) to provide a good 
degree of freedom upon which to organize the perturbative spectrum of K to any order (as 
usually happens in degenerate perturbation theory after degeneracies are lifted at one-loop). 

For our immediate purposes, eq. (3.19) constrains two-loop color structures. One easily 
sees that no double-real color structure satishes it by itself: for example, using (2.18), the 
hrst line of eq. (3.11) linearizes to give terms which replace three Reggeons WtWjWk by a 
single one VPq. Cancelling this term uniquely hxes the range-three part of the single-real 
contribution (to the form in eq. (A.10)). From other terms one constrains the double-virtual 
and range-two kernels. In this way, using in addition that double-real terms are signature- 
even, we hnd that the two-loop Hamiltonian can be parametrized by at most three angular 
functions: 


iL(2) = 


^ fd^d^ i2) abc( {LfUS'--Rf){Lp^,?-R^.)Rl 
dvr dTT l^-(L“-t/o-'Rf)(L5-<RpL- 




d d ( J^abc pa'h'd jrbh'jjcc' ^(Tjaa' i jjaa'\ \ / t a r^a' , r^a' t a\ 




(3.20) 




(2) (2) 

The hrst one is shown in hg. 7. Since ATV/qq, have already been determined 

from double-real emissions, effectively eq. (3.19) predicts all virtual corrections, up to a term 
proportional to the leading-order structure (the last line). The physical interpretation is that 
gluon Reggeization entails nontrivial real-virtual connections, which was indeed the original 
observation [11, 12]. 

In appendix A the prediction (3.20) is compared with the direct evaluation of single- 
real terms eqs. (3.15), (3.18). It turns out that there is a subtle loophole in the above 
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Figure 7. A combination of real and virtual color structures allowed by Reggeization. 


argument: in the non-global log context, U°‘^ —>• is not symmetry but only need to send 

the Hamiltonian to its complex conjugate. Thus signature is not conserved. The ansatz fails, 
by a single signature-odd term: 



OLQjaQk 


O^ik 


log 


aoiOofe/ oiQiaoj 


LfLpS'' Rl +Ll t/n" 


” TiCL 
itl itj 


(3.21) 

The origin of this term is simple: the imaginary part of the one-loop soft current (A.l). Its 
physical significance will be discussed shortly. The explicit computation in appendix also 
yields the yet-undetermined signature-even function: 


^( 2 ) _ 


aoiaoj 


( 2 ) I L 1 “b 

-yy -b bo log ^ 
daojcrojj 


( 2 ) MCa 


20ni?Tp 

9 


8 nsTs 

9 


(3.22) 

This contains only the two-loop cusp anomalous dimension [34] (normalized to 7 ^^^ = 1) and 
the one-loop /3-function.^ Physically, the bo term is fully dictated by the collinear anomaly 
discussed subsection 3.5, while the cusp anomalous dimension term can be checked to provide 
the correct Sudakov double logarithms in the limit of a narrow jet cone. 

How to explain the real-virtual connections (3.20) from the perspective of non-global 
logarithms? Perhaps one could use the Feynman tree theorem [36] . This is a way of computing 
loops by putting one (or more) propagator per loop on-shell. Indeed, putting a gluon on-shell 
in fig. 5 one can recognize diagrams of fig. 4, so at least schematically this seems to work. The 
tree theorem was streamlined and generalized to higher loops, with at least partial success, 
in refs. [37, 38]; it would be interesting to see its implications in detail. 


3.4 Double-virtual terms 

Let us now make sure that the ansatz (3.20) is not missing any virtual corrections. A priori 
these could involve two color structures: 

ir'" E {hrjkKR)Rl - KjkR^R^jRl) + E + KjLiq). (3.23) 

i,j,k i,j 

The coefficients are constrained by the KLN theorem (vanishing for = <5“^) and by Lorentz 
invariance, but without signature, unfortunately this does not fix them uniquely. (Below we 

® The appendix uses the so-called dimensional reduction scheme. In conventional dimensional regularization 
(CDR), more commonly used in QCD, a simple coupling redefinition [35] gives: —>■ . 
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give an example of a signature-odd function satisfying both constraints.) We resort to explicit 
computation. 

The two-loop soft anomalous dimension is known to take the ‘dipole’ form [5, 39, 40] 


r(2) 
^ n 






-2pi-pj - iO 


+ collinear terms. 


(3.24) 


This gives the divergence of the amplitude after subtracting the square of Since we are 
instead subtracting we need again to switch to the collinear-subtracted barred scheme 
(A.5): 

r^2) = r^ 2 ) ^ + collinear terms . (3.25) 

We omit ‘collinear terms’ which depend on only one leg at a time, since these are trivial to fix 
using the KLN theorem. Concentrating on the terms which have nontrivial color structures 
and which are not so easily fixed, the calculation of eq. (3.25) is rather straightforward and 
detailed in appendix A. The outcome confirms that no additional terms besides (3.21) need 
to be added to the Ansatz (3.20). 

We can now interpret this term. First we observe that it can be mostly removed by a 
finite scheme transformation. Namely, if we set 

^ g ^ log(aij)(L“L“ - + 0{al, c) , (3.26) 


where the MS density matrix is the minimally subtracted one we have been working with 
so far, then the two-loop Hamiltonian in MS^ gets shifted by a commutator with and a 
/3-function term. It is easy to check that the commutator term precisely cancels (3.21). The 
/3-function term then replaces it by 

(3.21) ^ 2mbo Y, log(«ii) - RiRj) ■ (3-27) 

This combination is Lorentz-invariant in an interesting way: under a rescaling of /3j, the j- 
sum become telescopic and simplifies to {Ci — Ci) = 0. This also satisfies the KLN theorem, 
being zero when L = R. The existence of this structure is the only reason we needed to use 
the explicit formula (3.24) to get the virtual corrections, otherwise the KLN theorem and 
Lorentz invariance would have sufficed. It violates the triangular structure (3.19) but since it 
is proportional to the /3-function this does not contradict the BFKL-based argument leading 
to it. It also has a simple and suggestive physical interpretation: effectively it replaces the 
spacelike couplings in the one-loop evolution, by timelike counterparts: 


5(/r2)i?“ ^ 5( V - i0)i2“, gip^)L<^^g{-p^ + iO)L\ (3.28) 


With hindsight, had we used timelike couplings in the one-loop evolution, we would never 
have had to write down eqs. (3.21), (3.26) nor (3.27). We will nonetheless continue to use the 
(more conventional) spacelike coupling. 
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3.5 Lorentz invariance and (lack of) collinear anomaly 


We have assembled all ingredients of the kernel, but we notice that the angular functions are 
not Lorentz-covariant: the arguments of the logarithms (3.12) are not homogeneous in /Iq, /3o' 
(and thus depend on the frame choice implicit in the normalization 13^ = (l,i;j)). This may 
be surprising given that dimensional regularization preserves Lorentz invariance. 

The simple explanation is that we did not write the one-loop evolution in a H-dimensional 
covariant form. What would constitute a Lorentz-invariant version is instead: 


/ ajj \ 
47r(27r)“2'^ yaoiaoj / 


47r(27r)“2^ 



+ 2e<5. 


( 1 ) 

ij;0 



(3.29) 


which differs at order e by the amount djj.g = 2 aQ "ao . The integrand is homoge¬ 

neous in all of jSi, f3j, /3 q and one may check that under a Lorentz transformation the Jacobian 
factor precisely cancels the change in the parenthesis. (The factor 4^ is for future convenience.) 

An 0(e) shift to an anomalous dimension, as usual, is equivalent to a finite renormaliza¬ 
tion of e.g. a scheme transformation. The density matrix in the ‘Lorentz’ scheme (3.29) 
is related to the MS one used so far, or better the MS^ scheme just defined, as: 


„ren,e ^ren,MS^ \ ^ 

a ' = a — / 

dvr 

(3.30) 

This shifts by a commutator as well as a /^-function term. 

This transformation is only well-defined because it contains both real and virtual terms: 
the middle integral in eq. (3.29) would otherwise be un-regulated even for e 7 ^ 0. This 
clash between Lorentz covariance and collinear divergences reflects the (now called) collinear 
anomaly of refs. [39, 40]. Here, the anomaly cancels between real and virtual terms. 

To achieve manifest Lorentz covariance we must still manipulate the expressions by using 
color conservation to add terms independent of some of the /?*, being careful with commutators 
as below eq. (3.3). Collecting these commutators is tedious but fortunately the task can be 
easily automated on a computer. We hnd (as it should) that the color structures in eq. (3.20) 
are preserved under these operations (see also ref. [41]). Only the angular functions change: 


d Op ( 1 ) 
dvr 


jraa / T a dcl i r>(i t cl\ Tara ryc 
^0 \^i^j +-^2 ^j) ~ ~ 


jya 

Kj 


a\ ren,MS 


(J 


r.{2)i _ ^(2)MS 


K. 


( 2)1 

ij-,00' 


= K. 


,, -I- Z y ^ ikflOjk - O ' + n . )j 

THjj-pp/ -|- 2A]i/j.pp/ 2Fji^-QiQ 

+ +FirM‘+Fjfm 


ij-,0 


(2)MS 

ij-,00' 

-m 

ij-,0 


^(2)lorentz _ ^(2)MS e(l) , 

“h ^00%;P 2 / dvT 


0 ' 


-FFij.Qff + F,j.q’„ + Fji.Q,„ — (0 o 0^)J-. (3.31) 

The functions Eij-^ and Tjj;p are arbitrary, with Eij-QQ/ = —Eji-fi/Q. The formula (3.33) below 

arises for Aj^ppi = - ^ -log and Eik-QQi = ^—— -log _ (With these 

choices the integrand on the last line vanishes.) 
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3.6 Final result for the evolution equation 

We record our final result for the two-loop Hamiltonian in the ‘Lorentz’ scheme (superscript £), 
which combines eqs. (3.20)-(3.22) with the hnite renormalizations (3.26) and (3.30). For con¬ 
venience we repeat the color structures, switching to the integro-differential notation (2.17): 

^( 2 )^ f 






T a TO Tjc jra jjb t c 


-t 





-t 


dvr dvr "bfc ;00 

d rig ^ ^0' j^{2)N=i,l ( fabc fa'b'c' jThh' ttoc' /rroa' , rTO,a'\\ / t o. jya' i Da'r i 

y j ^o' - +^o' 

C(ij 




dvr aoiaoj 


{nyL'} + L^Rj) + (3.32) 

Here = ^ = ^ {Lf'US'<^ - Rf), = {Ug-'Rf - Lf) and = f 

the color rotations L and R being differential operators dehned in eq. (2.17). All products 
of L“’s and 7?f’s are implicitly symmetrized and normal-ordered to the right of Uq, Uq/. The 
third term is simply the one-loop result (2.Id) times the cusp anomalous dimension (3.22). 
The angular functions are: 

^(2)e _ oii- ^ 

ao*ao'iA:^jfc;00' ~ 


Hj OiO'iOiO'jC^Qk _j_ ^ikOijk CKifcCtQ'j^ OiQ'iOijk O^jfcQ^OiQ^OO'Q^O'fc 


aoo' 


aoiaoja^if^ aokOio'k cnjkaoiao/k aoo'ao'fc 


OiQk^O'iOiO'j 


OiikOiQj . ajfcaoqaoo'Oofc , OiikOio'j . aifcaoo' 
log-^ log 


K. 


(2)7V=4,£ 


ijfiO' 


aofcaoo' 


aoi«oo'«o'i 


O^Q/k^Oi^^Oj 
aijaoo' 

2 log — -h 

OiQ'iOiOj 


ao'feOoo' oiokOio'i 

aijaoof 


O^OiC^jk OfjfcOiQO' 


1 + 


aoiO^o'j — aQuaoj 


aofcaoo' oio'kaoj 


log 


OlQ'iOlOj 


(3.33) 


This is the complete result in A/" = d SYM. In a general gauge theory with np flavors of Dirac 
fermions and ns complex scalars in the representation R, there additional contributions from 
matter loops, also obtained in eq. (3.12). Upon restoring group theory factors corresponding 
to representation R, in accordance with the square of fig. d(b), these can be written: 


j^{2)N^A ^ 



'1,3 


d^Q dClQi 1 
dvr dTT CKoo' 


aij log 




aoi«o'i “ o'o'iO'Oj 






X [inpTrRlT^UoT'^’ul,] - uf U^' - {urTr - 2Ca){U^^')} 

{LfRf+RfL-) 


f f dHo dUo' 1 

aoiOo'i + ao'i«Oi , aoicxo'j „ 

Ji,jJ 47r dvr 2a2g, 

io^ ^ 

aoiOoij — aoHOioj ao'i«Oi 


+ 


I 2{ns - 2nir)TrR[T“C/or“'c/J,] + 2r^^r'^'^'U^' 
\ -{{ns- 2nF)TR + CaW ^) 

2mbQ log(ay) (L“L“ - R^R^^). 


(3.3d) 




All sums are individually Lorentz-invariant (invariant under rescalings of the individual /3j). 
The first term is the contribution of two chiral N = 1 multiplets (minus the four adjoints in 
Af = A SYM) and the second term collects remaining scalars; bo = ^(IIC^ — AurTr — usTr). 
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4 Comparison with BFKL and conformal transformation 


As mentioned in the introduction, the same Hamiltonian K governs the Regge limit. Hence 
the reader familiar with the literature on the Regge limit, in particular the Balitsky-JIMWLK 
equation, will have recognized several equations by this point. Let us now discuss the con¬ 
nection in detail. 

Physically, as sketched in the introduction, the connection originates from the existence 
of a conformal transformation which interchanges the x"*" = 0 light-sheet and future (null) 
infinity. This interchanges the target residing at x"*" = 0 with the color rotations in the 
definition (1.2) of cr[U]. It is given explicitly as [17, 42, 43] 


y = 




y± = 


x± 

/iX+ 


X, 


y =x 


2x+ 


(4.1) 


where /x is a reference scale. This maps the Minkowski metric ds^ = —2dx~^dx +dx\_ to 
a multiple of itself, as one may verify. Points approaching the BFKL target, x+ —0, are 
mapped to infinity along the null direction oc (z?*^,/3 _l,/I z) = ( ^^^^^ > “^xx, ^ ^^^ )■ In 
this way the transverse plane of the BFKL problem is mapped stereographically onto the 
two-sphere at infinity of the non-global log problem. Provided that the conformal transfor¬ 
mation (4.1) preserves the Lagrangian, this map predicts that K should go into the BFKL 
Hamiltonian upon substituting [17]®: 


(Xij 



Xj)l , 



d‘^x± 

TT 


(4.2) 


We now verify this equivalence directly, beginning with the case of AA = 4 SYM where 
conformal symmetry is unbroken. In the general case we will find a discrepancy proportional 
to the /3-function, as anticipated. 


4.1 Comparison in AA = 4 SYM 

It is instructive to consider a special case: we act with on a dipole U 12 = Tic[UiUl]. The 
form (3.33) is particularly convenient for this since vanishes when i = k ov j = k. 

The only terms in the first line are thus Furthermore the remaining 

lines simplify since = 0, and one can check that 


R'lSp _ _ _ j^i2)N=4,£ _ ^(2)7V=4,£ 

'''■ 112 ; 00 ' ~ -^'■ 221 ; 00 ' ~ -^'- 12 ; 00 ' -^'■ 12 ; 0'0 


(4.3) 


In this way all two-loop color structures in the dipole case are expressed in terms of a single 
angular function. To evaluate the color factors we recall that while LfUi = in the 

® Here we use a normalization /3° A 1 which differs from that adopted elsewhere in the present paper and 
in ref. [17]. This has not effect in Lorentz-covariant expressions such as eqs. (3.33)-(3.34), but one should 
remember to include additional factors /3° = (1 -f if using non-covariant formulas such as eq. (3.12). 
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antifundamental one has that L 2 U 2 = —U^T^' (this ensures that (Li+L 2 )C/i 2 = 0). Writing 
if^-bcj'c _ coiiecting terms one easiiy finds 


ii'( 2 )Tr[C/iC/j] = 2 


(PQq jA2)N=4:,e 

dvr dvr 12;00' 


rxTjciCL^ jjbb' jjclcl' jjby jjcicl' jjby 
ZUq UQf — Uq Uq — Uq, Uq, 


X (^Tr[[T“,r^][/ir“'T'''c/]] +Tr[r^r“C/i[r^',T“']f7]]) 


+ 


dvr^Cyi f d?Qo ai2 


dvr aoiao 2 


Tr [r“' UiT^ul] i7o“'“ - CaTi') . (4.4) 


This formuia is identicai to the conformai form of the two-ioop evoiution obtained by Baiitsky 
and Chiriiii, eq. (6) of ref. [44], with l^ere ~ “i^lthere’ expected.^ In the 

pianar iimit eq. (4.4) can be reduced to a ciosed noniinear equation for a dipoie function; we 
refer to the iiterature for further detaiis. 

Rapidity evoiution for generai products of Wiison iines in the Baiitsky-JIMWLK frame¬ 
work has been obtained recentiy [41, 45-47], extending eariier resuits for two [44, 48] and 
three Wiison iines [49, 50]. Given the mutuai agreement between these works, here we oniy 
compare directiy against the conformai form of ref. [41]. Since the stereographic projection 
identifies the SL{2, C) conformal symmetry of the transverse plane with Lorentz symmetry 
of the two-sphere, this should match with the Lorentz scheme here. 

The comparison is in fact straightforward: the range-three kernel /C 3^2 shown in eq. (5.12) 
of ref. [41] is literally the first four terms of our remaining two terms in 

arise from the telescopic term F in eq. (3.31) hence do not affect the range-three part. (These 
terms are helpful to manifest the convergence at /3o —>■ /3o'-) Furthermore, the integral repre¬ 
sentations for /Csq and JC^fi in ref. [41] reproduce the real-virtual pattern embodied in the first 
line of eq. (3.32). This demonstrates the agreement of range-three interactions. Combined 
with the agreement in the dipole case, this establishes the complete equivalence of eq. (3.32) 
with ref. [41] (and thus, by extension, refs. [45, 46, 49]). 

In principle, upon linearizing around U = 1, one also expects complete agreement with the 
interactions between Reggeized gluons obtained in the BFKL approach. For two reggeons the 
agreement was demonstrated at the level of eigenvalues [44, 48, 51, 52]. For three reggeons, 
it was noted in ref. [49] that a scheme transformation appeared to be lacking in order to 
match with ref. [53]. This issue should be clarified further. Here we simply note that there 
is a natural candidate: the next-to-leading order inner product (correlator of Wilson lines) 
[54, 55]. In the BFKL approach the inner product does not receive loop corrections (the 
transverse part of the Reggeon propagator remains l/p^), so only after this effect is removed 
by a scheme transformation, should agreement be expected. 


Since subtraction terms are not written in exactly the same, there is an apparent discrepancy: 


■I 


FQo j^(2) 


Itt Itt 


Zo' {ur'uZ - Tr[[r“,T'’][/i[T“',T*’']l/t] cx J ( 


[ _ j^(2) 


;00'J ■ 


However the integral vanishes, proving the agreement. This is easily shown by noting that being absolutely 
convergent, the integral defines a Lorentz-covariant function with the same homogeneity in /3o,/3i,/32 as the 
integrand, hence must a constant times q, ■ The constant vanishes by antisymmetry in (j3i A)- 
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It is interesting to compare technical aspects of the calculations. The tree-level soft 
current (3.1) is reminiscent of the light-cone gauge amplitudes in eq. (43) of ref [48]. The 
subtraction of sub divergences in eq. (3.10) is similar to the -|- prescription derived in refs. [48, 
56]. The transformation to the ‘Lorentz scheme’ (3.30) is identical to that leading to the 
‘conformal basis’ in refs. [41, 44]. As a significant technical simplification, however, we saved 
the Fourier transform step. Also the reliance on standard building blocks made it possible to 
benefit from results in the literature, namely the soft currents and collinear splitting functions. 


4.2 Comparison including running coupling 

Having demonstrated the agreement in A/" = 4 SYM, let us now compare the fermion and 
scalar loop contributions to the Balitsky-JIMWLK and non-global logarithm Hamiltonians, 
e.g. the terms proportional to np and ns in eq. (3.34). Performing the comparison with 
refs. [44, 57] we find that the two Hamiltonians agree for the most part, except for the 
following discrepancy (setting Zij = Zi — Zj): 


arj 



'1,3 


dPzo 


TT 


(j a T a X TDa pi 


-27Tibo / log(4.)(4“i“-^“^“) 


'1,3 


z 2 


log(/i^4) + ^^.2 . 2 log ^ 


z2 z2 


Z, 


OJ. 


(4.5) 


where as before // is the MS renormalization scale. In particular, the difference is proportional 
to the first /3-function coefficient, as predicted [17]. 

The origin of the discrepancy (4.5) is clear; the inversion y~^ —)• in (4.1), which 

relates the BFKL and non-global log Hamiltonians, is only an isometry up to the Weyl 
rescaling dSy —)• {yLy~^)~‘^dSy. This is not a symmetry in a non-conformal theory. Physically, 
BFKL and non-global logarithms describe infinitely fast and infinitely slow measurements of 
an object’s wavefunction, which would not normally be expected to be connected without 
conformal symmetry. 

For future reference, we note that a general theory deals with Weyl transformations in 
non-conformal theories (see for example [58]). The essential feature is that, starting from the 
BFKL side and performing the conformal transformation (4.1), one ends up with a coordinate- 
dependent coupling constant: 

5" = / ’ oisif^ony^) = asino) log(/xy+) + . (4.6) 

In other words, the BFKL Hamiltonian in QCD in principle controls non-global logs in QCD 
but in an imagined setup with a coordinate-dependent coupling. Contrary to real QCD, in 
this setup a narrow jet never hadronizes: the increasing coupling due to the growing size of 
a jet, is compensated by its falloff at large y~^. Thus effectively the coupling is set by the 
angular size. This reflects that angles map to distances in the BFKL problem. We will not 
pursue eq. (4.6) further here, but in any case it is clear that to all orders in perturbation 
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theory the difference between the BFKL and non-global Hamiltonians must be proportional 
to the /3-function (up to finite scheme transformations). 


5 Higher loops and exponentiation 


It is instructive to extend the general analysis of section 3 to higher loops. We will (mostly) 
ignore collinear subdivergences here, concentrating on the soft divergences. 

We can organize terms according to the number m of wide-angle partons {U matrices) 
added to an underlying n-jet event. Our starting point is the known exponentiation of virtual 
corrections (2.1), which gives the m = 0 case: 


(70[U] 


Pe-fo Tn = pg- fo^ Ko ^ren 


(5.1) 


For the next case of one wide-angle gluon, a formula was derived in eq. (3.14). We reproduce 
it here, in abbreviated notation, omitting U matrices, the angular integration, energy 

measure, and absolute value squared on the matrix elements: 


o-<i 


= : Ve- ^0 ^0 


<7n 


V) + 


Si{a-,a)Hnip) -f 


Hn+l (O) 1*) 


0<a<iM 


fi<a<oo 


(5.2) 


The colons instruct us to normal-orders terms according to their renormalization scale (largest 
argument to the right). As in subsection 3.3, the first integral is identified as a shift A"i(/r) = 
—S'i(/r;/x) to the exponent. The remaining (finite) terms then define the hard coefficient 
(T<{‘(/u). Moving on to two real emissions, we follow eq. (3.9) and write the cross-section as 
independent emissions plus an additional piece: 


a<2[U] = : pe“^o'‘(^o+i^i) 


^<T(l^) + 


53 ( 0 , 6 ; b)Hn{p) + 


Hn+2{a,b-,^) 


0<a<b<fi 


fi<a<b<oo 


+ 


{Hn+2{a, 6 ; /r) - 5i(a; a)Hn+i{b; p)) 


0<a</.t 

fi<b<cx 


(5.3) 


We have introduced the ‘connected’ squared soft current defined as (all factors being evaluated 
at the same renormalization scale) 


52^(a,6) = 52(a,6)-5i(a)5i(6), 

53 ( 0 , 6 , c) = 53 ( 0 ,6, c) - 5i(a)52(6, c) - 52 ( 0 ,6)5i(c) -f 5i(a)52(6)53(c), etc. 

(In the present abbreviated notation we recall that each factor is a squared soft amplitude, 
Si = |5ip.) Again the first integral in eq. (5.3) is identified as a shift to the exponent, 

= - [ 53 ( 0 ,/r;//), (5.5) 
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which generalizes eq. (3.10) to include virtual loop effects to all orders. The (finite) remainder 
then defines • 

Using this method it is straightforward to extend the calculation to more radiated par¬ 
ticles. For three radiated particles, for example, after pulling out 'Pe ~we find 
again that particles with energy > /r decouple from divergences, no subdivergences, and the 
following shift to the anomalous dimension: 

Kz{lj)= -j . (5.6) 

0<a<b<fi 

The absence of subdivergences (finiteness of as e —)■ 0) is manifest from the fact that 
the expression involves only connected squared amplitudes, which vanish near the boundaries 
a —)• 0 or a, 6 —7- 0. This itself is a consequence of factorization, or more precisely eq. (2.5) in 
the form 


lim 5 '(q,i, ..., Uri) A*) — *S'(ui, • • • j Ofcj h')'5(Q.fc-(-i, •.., dni A*) • 

It is now clear how to generalize the pattern to higher orders. In fact a simple guess appears 
to give the anomalous dimension K to all orders: 


(5.8) 


The exponential factor has a simple physical interpretation as an ‘exclusion time’ effect, and 
we recall that a’s are the energies of real radiated particles. We have verified explicitly (with 
the help of a computer) that exponentiating K using eq. (2.7) reproduces all contributions 
where up to at least 9 real particles have energy below /r, so we believe that the formula is 
correct to all orders. 

Equation (5.8) is one of the main results of this paper. It expresses, to all loop orders, the 
Hamiltonian governing non-global logarithms as a convergent integral over finite, well-defined 
building blocks, generalizing the eqs. (3.10) and (3.15) used in the two-loop computation. 
The building blocks are the squares of the infrared-renormalized soft currents (which include 
virtual loops to all orders), defined in eq. (2.5). Only the part of the renormalized currents 
are needed, in agreement with ref. [59]. 

Since the exponent K is manifestly finite as e —0 (being expressed in terms of connected 
squared soft currents), the formula also demonstrates to all loops that infrared divergences 
exponentiate according to eq. (2.7). The physical inputs were the known exponentiation (2.1) 
of virtual corrections, and the factorization of successive real emissions (5.7); eq. (2.7) is seen 
to be a purely combinatorial output. 

To fully prove eq. (2.7) one should address the issue of collinear subdivergences, omitted 
in the present discussion. Physically we expect these to cancel, since (j\U] is collinear-safe. 
In subsection 3.3 this was made manifest by defining collinear-subtracted real and virtual 


= r„ + 


oo 

E 

k=l 


-pgJai 


K 


• • • ) ^k—li A^j A*) • 


Oj\ <C. . "C/J- 
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contributions, such that their sum was unaffected by the subtraction. We have no reason to 
think that this couldn’t be achieved at higher orders as well, following the method of ref. [25]. 

We should mention that eq. (5.6) gives the evolution equation in a non-minimal scheme: 
the two-loop exponent in eq. (5.5) depends on e through a factor a~‘^^ and thus differs 
from the MS result of this paper by terms proportional to e. These are interpretable as a 
renormalization which shifts by a finite commutator, giving, instead of eq. (5.6): 

-J - I ^(Kj')52^(a,M) + 52^(a,/i)i^f^)log^ + 0(e). (5.9) 

0<a<b<fi 0<a<fj, 

(' 2 ') 

Discrepancies in at 0(e) add other commutator terms. All these are unrelated to a fur¬ 
ther finite renormalization needed to make Lorentz covariance manifest. As in eq. (3.29) it can 
be fully predicted by upgrading the two-loop result (3.33) to a D-dimensional covariant form. 
Although these finite renormalizations become combinatorially very complicated at higher 
loop orders, being finite they cannot interfere with the statement (2.7) of exponentiation. 

Finally, we list the ingredients which enter eq. (5.8) at three-loop: 

• The tree-level soft current for three soft gluons S^\ at one-loop for two gluons 

(2) (3) 

and two-loop for one gluon 5} , and the three-loop soft anomalous dimension Tii b 

• The next-to-leading order 1—7-2 and tree-level 1—7-3 collinear splitting amplitudes [60] . 

The two-loop soft current and three-loop soft anomalous dimension are presently known for 
two hard partons [61-63]. Unfortunately this will not suffice for non-global logarithms nor 
BFKL in general, since each radiated gluon counts like a hard one from the point of view 
of softer radiation.® However, for dipole evolution in the planar limit, everything needed is 
known. 

6 Conclusion 

In this paper we considered a ‘color density matrix’ which aims to characterize soft radiation 
in gauge theory. We argued that it should resum large logarithms arising in the presence of 
wide-angle phase space cutoffs, so-called non-global logarithms, to all orders in logarithms and 
l/Nc- We proved the all-order exponentiation of infrared divergences for this object in terms 
of an anomalous dimension K (see eq. (2.7)), constructed formally in eq. (5.8), modulo one 
technical assumption regarding collinear subdivergences. We explicitly computed K to two- 
loop (eqs. (3.32)-(3.34)) and performed a number of checks on this result. We also stressed 
the equality between K and the BFKL Hamiltonian, which allows our results to be viewed 
as an independent derivation of the next-to-leading order BFKL Hamiltonian, obtained here 
directly in a novel, compact form. 

® An interesting possibility is that other constraints, such as the KMS condition; Lorentz symmetry 
(SI/(2,C)); collinear singularities; CPT symmetry; and Hermiticity of the BFKL Hamiltonian could uniquely 
fix without these building blocks, effectively determining them. 
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The procedure to calculate a cross-section receiving non-global logarithms was sketched 
in the introduction. One distinguishes infrared and ultraviolet scales, which are to be con¬ 
nected by evolving using K. At both ends lie finite quantities: an TR measurement’ which 
contains details of the experimental definition of a ‘soft’ particle, and corresponding vetoes; 
an ‘UV measurement’ which depends on the initial state and possible vetoes imposing hard 
jets in the final state. The logic of factorization being that their calculations are independent 
of each other, we focused in this paper on the (universal) evolution K. Study of the (expect¬ 
edly) finite, but process-dependent, measurement functions is left to future work, as well as 
phenomenological studies. 

Mathematically, K is an integro-differential operator acting on functionals (t\U] of a two- 
dimensional field of unitary matrices U{6) (e.g. SU{3) matrices in QCD), with 6 an angle in 
the detector. This means that K cannot be diagonalized explicitly. Although it is a quite 
complicated object, it is a useful starting point for further approximations. These include, 
as reviewed in section 2, numerical Monte-Carlo techniques at finite Nc, reduction to an 
ordinary integro-differential equation at large Nc, or linearization d la BFKL around U = 1. 
We hope that one of the compact forms obtained in this paper will prove convenient for a 
next-to-leading order numerical implementation. 

For application to hadron colliders it will be important to go beyond the limitation of an 
initial color-singlet object, as done in this paper, and allow for initial state radiation. This 
could lead to additional (super-leading? [64, 65]) effects related to subtle color-dependent 
phases in collinear limits [66, 67]. 

The formalism does not distinguish between global and non-global logarithms, but it is 
easy to see how it simplifies in the case of global observables. For example, when radiation 
is excluded everywhere but inside narrow cones, the IR averaging procedure sets {U) = 0 
outside these cones which effectively shuts down the real terms in the evolution. It is then 
dominated by virtual effects, as is usual for global observables. It is only for observables 
sensitive to details of wide-angle radiation that the complications of the formalism kick in. It 
would be interesting to connect the present approach with that of ref. [68], which deals with 
recursive infrared and collinear safe event shapes (‘rlRC’). 

There has been recent activity regarding formal aspects of measurements at infinity, in 
connection for example with the Bondi, van der Burg, Metzner and Sachs (BMS) symmetry 
[69, 70]. The density matrix construction could be useful in this context. 

From a theoretical perspective, the Hamiltonian K connects, in a unified way, the fol¬ 
lowing gauge-theory concepts: the cusp anomalous dimension (governing global logarithms); 
the KLN theorem (cancelation of collinear and infrared divergences); the factorization of soft 
radiation; the BFKL equation. 

The equivalence with BFKL, verified explicitly in section 4, is a consequence of conformal 
symmetry [17] and is an equality up to /3-function terms (fixed by comparatively simpler mat¬ 
ter loops (4.5)). The basic physical intuition is summarized in fig. 2. Remarkably, properties 
manifest in one context are not necessarily so in the other. 

For example, one fundamental assumption in both the BFKL and Balitsky-JIMWLK 


- 28 - 


frameworks is that transverse integrals should be saturated by transverse scales ~ t <C s, 
ensuring that rapidity logarithms (logs) arise only from longitudinal integrations [18, 71]. 
While reasonable it is unclear how one would prove this directly beyond the current state of 
the art, e.g. next-to-leading log. The correspondence with non-global logarithms immediately 
implies it to all orders, since it amounts to the amply understood cancelation of collinear di¬ 
vergences. The non-global logarithm formulation may also be computationally advantageous, 
as discussed in sections 4.1 and 5. 

In the other direction, the phenomenon of gluon Reggeization suggested a compact way 
to write the evolution equation (see eq. (3.20)), which manifests a connection between real 
and virtual effects. Intriguingly, we found that these relations could perhaps also be explained 
by the Feynman tree theorem. It would be very interesting to see if either of these approaches 
generalizes to higher loop orders. 

Finally, we mention that the simplest non-global logarithms to resum in this frame¬ 
work (beyond the planar limit) involve situations close to the linear regime 17 ~ 1, where 
the linearized equation has lowest eigenvalue the well-known Pomeron intercept — 

Naively this regime might correspond to multiplicity-type measurements, e.g. counting away 
jet charged tracks as a function of angle and an energy cutoff. Perhaps this or some other 
observables could provide an indirect experimental handle on the BFKL Pomeron. 
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A Single-real and donble-virtnal contribntions 

We detail our evaluation of eq. (3.15). The starting point is is the one-loop soft factor, which 
we take from refs. [60, 72, 73]: 

5«'^“(A:i) =g^Yl + ^9^ E (A.l) 

i i,j 


with 

sP^ik) = 


6 /3i-k ’ 



P i-2(3,-k-i0)i-2(3yk-i0) \ 
2y/3j-k /3i-kJ V p{-2/3i-l3j - iO) J 


This is transverse and consistent with eq. (2.6). We need to convert this to the barred scheme, 
SP = where , (A.2) 


- 29 ~ 









where H implements the subtraction in eq. (3.17) of collinear splitting functions. Since the 
splitting functions for all but the radiated gluon cancel in the commutator, we will only need 
the gluon splitting function 


|Splitg(a/3o,6^o')|^ 


CA{b - a)-^^ 
ab{b~‘^^) 


x(l—x)aoo' 


+ - 


(^Weyl 

«00' 


fa;(l-x)(nfJ-2n^Jyi+2)/ 


(A.3) 

where x = a/b. The prefactor has a kinematical origin and accounts for the change in the 
measure b^~‘^^db. The computation of such functions is standard [32]. In the x-dependence 
one can recognize various DGLAP kernels )(x), as expected. We use the dimensional 

reduction scheme so the parenthesis does not depend on e. (Regarding color factors we 
recall that we show intermediate formulas only in a theory with color-adjoint matter.) The 
scalar contribution to the splitting function is polarization-dependent and for us the most 
useful information will be its dot product against divided by Pi-Pj'. this is what enters 

eq. (3.15). This is given by 


f{Po,Po']Pi,Pj) 


[aoi — ao'j] [aoj — ao>j] 

®b'®oo' 

[ao>i{aoj — aoo>) — ao'j{aoi — aQo')]"^ x ■ 

- 2 -1" convergent or telescopic . 

‘2o'ijCXQQiO'(j'iO'Q'j 


(A.4) 


The first form is obtained directly from the Feynman rules and makes manifest that the 
dependence on Pi,Pj is consistent with factorization. We will prefer the second form, which 
provides a closer match with eq. (3.12) and also yields a simpler integrated expression. It 
differs by terms which are either convergent or vanish using color conservation. Computing 
the integral in (3.17) we then obtain 




RiRj{L 2 {aij)+ log4loga^) 

ij i=g ^ 


+ 7^^ +6olog^+^log2 



(A.5) 

with L 2 {x) = Li 2 (l — x) ^ log^(x) — The sum runs over gluons to stress that we haven’t 
computed the other cases, and the cusp anomalous dimension is in eq. (3.22). 

The commutator then easily yields 




+ {bo - 27riCA) log 






(A.6a) 


Pj-Po 


2 ( 



+T2(c7j)—T2(aoi) — -^2(aoj) 


otQiaoj \ 


(A.6b) 


We stress that only the O(e^) terms of were needed to obtain this. It is noteworthy that 
the —tt'^Ca/Q from the original soft function, the —tt'^CaI^ from the scheme transformation. 


- 30 
















and the —Ca log(e“*’^)^/2 from the phase of the logarithm have nicely canceled to leave the 
cusp anomalous dimension. 

Substituting into eq. (3.15) the soft factor produces two color structures 




linear in U 


E 

* ' 


- L]Llur'R^ 

+ E / (i? < + . ( A.7) 


which multiply the angular functions 

G'lh = -4Sf , Gg = . (A.8) 

These can be evaluated explicitly using (A.6). The remaining linear-in-f7 contribution, the 
subtraction (3.18), is simply 


= — 

'^sub. 
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d^Qo d'^flg' OCijCA 


dvr dvr 




(T“i?“ +Rf m 


-d 

''■Weyl ^ 


“00' 


+h"f'-2"wll+2)/ 


(A,9) 

This is to be added at the integrand level to eq. (3.12) and removes its collinear divergences. 
(Since the cancellation occurs at the integrand level in eq. (3.10), it is justified to set e = 0 
in both.) 

The result (A.8)-(A.9) is now to be compared against the prediction from the ansatz 
(3.20), which gives the same color structures and predicts the first angular functions as 


^(2) I 

2J'7c; 0 I ansatz 
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dTT 
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( 2 ) 


ikj]00' 
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\aoiaok aoiOioj 


)G 


5 log 


V ^jk ) 


+ L 2 {ajk) — L 2 {aoj) — ^ 2 ( 09 ^)^ . 


(A.IO) 

This agrees precisely with eq. (A.6a), up to the iir term recorded in eq. (3.21). For the other 
structure 


^( 2 ) I _ I 

; 0 ' I actual I ansatz 
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OlOiOiQj 
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4 ^ ^ij-,00' 
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daoiCtoj 


- A 




(A.ll) 


( 2 ) 

which fixes Ab.g, as recorded in the main text. 

Finally we check the double-virtual terms. To get the prediction from the ansatz (3.20) 

we need to integrate (A.IO). The T 2 (aofc) terms look scary, but they cancel out trivially 

( 2 ) 

because one needs only the total antisymmetrization of G)jfc.oo' modulo terms which do not 
depend on all three labels simultaneously. The integral is still a bit nontrivial but we could 
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simplify its antisymmetric part using integration-by-parts. We omit the details and quote 

( 2 ) 

only the rather simple result for the contribution, 

ir'^RIR'-RiY, j ^ (A.12) 

iijik ^ '' 

Finally the other term in the ansatz is (dropping terms depending on one particle at a time) 
- (iP log(aij) - bo{L2{aij) + log4:logaij)) . 

The preceding two equations are easily verified to be in perfect agreement with the commu¬ 
tator (3.25), proving that the ansatz does not miss any double-virtual term. 

As a final comment, we note that the L 2 function and most log2’s have a simple origin: 
the scheme change (3.30). For example / if^ 2 w)-°cr ^ijfi = -^ + '^+L2{aij)+\ogA\og{2aij) + 
0(e). With hindsight, we could have saved ourselves much algebra by switching from MS 
to the Lorentz-covariant scheme in the first step, which would have prevented L 2 from ever 
appearing. 
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